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I. INTRODUCTION 



In recent years it has been shown that one can create a condensed matter/gauge theory duahty using the AdS/CFT 
correspondence This duahty can give us insight in how a strongly coupled system behaves by reformulating 
the theory in terms of a weakly coupled theory of gravity. As of now the Hall effect Nernst effect and 
superfluidity /superconductivity H, [1| have been explored. Some reviews are [l-Q- 

Most of the work has been done with superconductivity. Initially the probe limit was explored, in which the 
chemical potential y^t — >■ whereas the charge of the scalar hair g — oo so that the product q/j, remains finite [sj . This 
work was later extended to include the full backreaction on the geometry of the system . Magnetic fields have also 
been added and found to exhibit the behavior of Type-II superconductors (lo| . Other gravitational duals have been 
able to describe systems exhibiting striped superconductivity (see, e.g., [HI). All dual systems show hallmarks of the 
strongly coupled, high Tc superconductors. 

In this paper we explore the properties of the system at the opposite end of the probe limit at which the charge of 
the scalar hair q is close to zero and the critical temperature Tc is small (small Tc/ ii limit) corresponding to a near 
extremal black hole. In Hartnoll, et al., discovered that a phase transition is allowed even for g = (see also 
discussion in [HI, [T^). It is not clear what the physical mechanism responsible for the condensation is near g = 0, and 
whether it is related to the one in the probe limit {q — > oo). To address these questions, we explore the near extremal 
limit by performing analytic calculations of the properties of the system. There is a critical line in the two-dimensional 
parameter space (A, g^) defined by [l^ 



2 3 + 2A(A-3) 



4 ' W 

which is determined by the AdS2 geometry of the horizon at extrcmality. Instability is expected above this line, with 
the critical temperature approaching zero as one approaches the critical line. By using trial functions, we show that 
indeed, there is an instability above the critical line ([T]), but the critical temperature vanishes on the line only for 
A > Ag. We compute the scaling dimension Ag analytically (numerically, Ag sa 0.64). For A < Ag, the critical 
temperature does not vanish on the line ([1]) and there is a phase transition below this line. We find an analytic 
expression for the true critical line (where Tc = 0) for A < Ag. Interestingly, for q = 0, we obtain instability down to 
the minimum value A = 0.5 determined by unitarity (at A ~ 0.5, we obtain Tc = 0). 
We also calculate analytically the energy gap 



Tr. 



E,^'-^^, (2) 



which is found to diverge as we approach the critical line {Tc — i> 0). A calculation of the conductivity (both analytic 
and numerical) shows that the profile of the conductivity exhibits a universal behavior near criticality. 

In detail, in section|lI]wc set up the system and fix the notation. In section lllll we determine the critical temperature 
analytically. In section IIVI we study the system below the critical temperature and calculate the energy gap both 
analytically and numerically. In section |V] we calculate the conductivity analytically near criticality and compare with 
numerical results. Finally, in section IVll we summarize our conclusions. 



II. EQUATIONS OF MOTION 



We are interested in the dynamics of a scalar field of mass m and electric charge q coupled to a J7(l) vector potential 
in the dynamical backgound of a 3 + I— dimensional AdS black hole. The action is 



where = dp, — iqA^. For simplicity, we shall work with units in which L = 1 and 16nG = 1. 
It is convenient to write the complex scalar field in terms of real scalar fields and 6 as 



Then the action reads 



R + 6- -F^'^'F^,^ - dp-^d^"-^ - q^-^'^id^e - A^f - mH'^ 



(3) 

(4) 
(5) 
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The phase 6 plays the role of a Stiickelberg field giving mass the the vector potential when the real scalar ■0 condenses. 
The local U{1) transformation reads 



(6) 



We can fix the gauge by setting 9 — 0. The form of the action ([5]) is more general than ([3]), because in the former, we 
can continue to negative values. This leads to a negative squared mass for the vector potential which signals an 
instability in flat space. However, in AdS space negative squared masses are allowed as long as they are above the BF 
bound. We shall allow to take negative values and find that there is a lower bound for a phase transition to occur. 
To solve the equations of motion, consider a black hole with planar horizon and metric 



1 



dz^ 
9{z) 



(7) 



We have chosen the coordinate z so that the horizon is at z = 1 = 0). The system is kept under fixed charge 

density g. We shall denote by p the value of the charge density in units in which the radius of the horizon is 1. 

All physical quantities are measured in units of the horizon radius, which hides an arbitrary scale in the system. 
This is remedied by reporting on scale-invariant (dimensionless) quantities only, as we shall be doing. 

Assuming that the scalar field is of the form ^(z) and the potential is an electrostatic scalar potential Aq = $(2;), 
the equations of motion are 



Y 


x' 


2 




.9 


2 


z 





2 



4 z 



-X 
3(.9 - 1) 



2 
2z2 



z'^g 

z'^g 
2 ^ 



* = 0, 



$ 0, 



0, 



= 0, 



where prime denotes differentiation with respect to z and to be solved in the interval (0, 1), where z 
and z = is the boundary. 

Near the boundary (z —> 0), we have g — 1, x — > and so approximately 



4' w * 



pz, 



where 



Ah 



± 



(8) 

1 is the horizon 
(9) 
(10) 



For rn^ < — | (Breitenlohner-Freedman (BF) bound [U), A± have an imaginary part and the system is unstable. 

For — T < < — 1> both modes are normalizable. While a linear combination of these modes is allowed by the 

.4 — . ^ . ... I — I, . . . 

equations of motion, it turns out that any such combination is unstable |15| . However, if the horizon has negative 

curvature, such linear combinations lead to stable configurations in certain cases [l^. Thus, the system is labeled 

uniquely by the dimension A = A± . For > — - 

that A = A_ > i (unitarity bound). 

Demanding at the horizon 



only the mode of scaling dimension A_|- is normalizable. It follows 



$(1)=0, 



(11) 



p is interpreted as the chemical potential of the dual theory on the boundary, p is the charge density on the boundary, 
and the leading coefficient in the expansion of the scalar yields vacuum expectation values of operators of dimension 
A±, 



The Hawking temperature is 



T 

71 



9'W 
47rv/-3''(0) ' 



= -x(l)/2 



(12) 



(13) 
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The equations of motion admit non-vanishing solutions for the scalar below a critical temperature Tc where these 
operators condense. 

Above the critical temperature (T > Tc), we have = and x = 0. The equations of motion reduce to 



= 0, 



£!$'2 _g^^ 3(.g - 1) 



= 0. 



(14) 



which yield the Reissner-Nordstrom black hole 



g{z) = l-[l + 



whose Hawking temperature is 



<J>(z) = p(l - z), 



^ 47r^ V 4 



(15) 



(16) 



III. THE CRITICAL TEMPERATURE 



To find the critical temperature Tc, we note that ai T = Tc the equation of motion for the scalar field ^ decouples 
from the rest of the equations of motion, because ^' can be treated as a perturbation near the phase transition, and 
consequently, g and $ are simply given by (jl5|) . Therefore, to find Tc, wc need to solve the equation of motion for ^ , 



Y 


2 


1'' + 




.9 


z 


52 ^2g 



* = 



(17) 



together with the boundary condition ^ z^ aX the boundary {z = 0) in the Rcissncr-Nordstrom black hole 
background ([15]). We will do so by working very close to the zero temperature limit of the black hole. At extremality 
we have 



. 



Near the horizon [z 



1), the scalar field behaves as 

[l-zf^ 



5± = A^^ 

2 V3 



3 + 2A(A -3) 



(18) 



(19) 



This suggests that there is a potential instability above a critical value of the electromagnetic coupling constant, given 
by the critical line ([1]) in the two-dimensional parameter space {/S.,q^), since both scaling dimensions 5± are complex 
and therefore no solution is regular at the horizon. We shall show that this instability does indeed occur and calculate 
Tc analytically. 

This instability has a geometric origin. Near the horizon the geometry is AdS2 x -v^rjth radius L2 = l/'v/O, and 
the scalar has an effective mass given by 



(20) 



The system is unstable below the BF bound for AdS2 space, 'rnl.^L\ < —\, which is in accord with ([T|). 

One expects Tc — ?> as we approach the critical line ([T} from above, whereas below the line no instability is expected. 
We shall show that this only partly true. For A > Aq, where Ag ~ 0.64 (we also compute Aq analytically), we obtain 
Tc = along the line ([T]). However, for A < Aq, Tc 7^ along this line and wc obtain instability below this line down 
to a true critical line (on which Tc — 0) that we compute analytically. Among other things, this implies that for g = 
we have instability down to the unitarity bound A = 0.5. 



First, we consider the special case m? ~ —2, in which A = 1,2, and a solution can be built on hypcrgcomctric 
functions. A similar approach for other values of m? yields Heun functions which are not as manageable. For the 
more general case, we shall rely on trial functions. 
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The critical temperature corresponds to a black hole near extremality, so we set 



^ = 3-e, I^ = To = -, (21) 



where e ^ 1. The critical temperature is given by 



Tr. Tf 







(22) 



To find the critical temperature (i.e., e), we need to solve the wave equation ()17p near extremality. It is a second order 
equation for a real function, therefore the general solution is given in terms of two arbitrary real integration constants. 
There are two boundary conditions to be satisfied, one at z = 0, where 4* ~ (A = 1, 2), and regularity at the 
horizon (z = 1). These two fix one of the integration constants, as well as the eigenvalue e. The second integration 
constant (an overall normalization) remains arbitrary, because the wave equation is linear. 

To solve the wave equation, split the interval [0, 1] into two overlapping regions, one away from the horizon, defined 
by 1 — z > e, and one near the horizon, defined by 1 — z < ^/e. These two regions overlap for e < 1 — z < ^Je, if e < 1. 

Away from the horizon (1 — z > e), we may approximate the metric by the extremal metric. Then we obtain the 
solution at the critical temperature away from the horizon. 



z - \z - zq J \z- zqJ y V3 V3 z- Zq J 

where zo,Zq are the two complex roots of g{z) with positive and negative imaginary part, respectively, and S± are 
the scaling dimensions p9|) with m? = —2. The multiplicative constant C is complex. We shall use the boundary 
conditions to determine its phase, whereas \C\ will remain arbitrary. 

Near the horizon (1 — z < \/e), we need to exercise care. Let us perform the coordinate transformation 

^ = 1 - ^C, (24) 

so that C < l/\/e- Notice that C can be large, and the horizon is at C = 0. 
Expanding around the horizon, 

9{z) = -^C5'(l) + I (-^C)' .9"(1) + • • • = jCil + C) + . . . , (25) 
where we used g'{l) = —AttT, g"{l) ~ 6 — 1 ) , and eq. (PT|) . the metric becomes 



ds^ = - 
6 



dC 



2 



C(i + C) 



df + ..., (26) 



where we omitted higher-order terms in e. One can verify that the Hawking temperature is still given by (|22l) . The 
electrostatic potential reads 

$ = ?C+..-, (27) 
and expanding in e, the wave equation for the scalar field at the critical temperature near the horizon becomes 

.2 c 



C(i + C)*" + (2C + i)*' + ^ 



1 + 



i + C 



^ = , (28) 



where prime denotes differentiation with respect to C. The acceptable solution at the critical temperature near the 
horizon is 

*ncar(C) = ^ (1 + Cr^/^'F i^-5+ ^, - ^; 1; -() , (29) 

given in terms of a different hypergeometric function, where yl is a real arbitrary normalization constant. The other 
(complex) normalization constant C will be related to A by matching with the solution ([25)1 far away from the 
horizon. Notice that the other solution to is discarded because it has a logarithmic singularity at the horizon. 
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It is easily deduced from the identity 

F[a, /3; 7; x) = (1 - xy-'--^F{^ - a, 7 - /3; 7; x) (30) 

that ^'ncar IS real, 

*„ca. = A{1 + C)^^/^ F (^-5^ + ^, -5+ + ^; 1; -^j = ^Lar • (31) 

Next, we match the two solutions in the overlap region. 

We approach this region from near the horizon for C 3> 1. The wave function becomes 



C*++c.c., (32) 



In the far region, we need to take the limit z — > 1. The corresponding wave function becomes 

2V2-i 



*fa,. ~ C(l - z^y- (^^— ^ ) (1 - zf+ + C.C.. (33) 

Using (|24p . we deduce the relation 

r(-<5_ + ^)r(-<5_-^) Ve^ Vi-^oV 



between the constants A, C, and e. Thus, the phase of C has been fixed by imposing one of the boundary conditions 
(regularity at the horizon). An overall normalization constant A remains arbitrary, as expected, since the wave 
equation is linear. 

The eigenvalue e will be fixed by imposing the remaining boundary condition at z = 0. Focusing on A — A-|_ = 2, 
we demand the asymptotics ^E* ^ near the boundary, giving 















) 











C(-zo)*- (^^j F\^-S_ + 2^^q,-S^-^;-26^:-2zoj+c.c.^0 . (35) 

By solving this constraint with C given in p4p (notice that the arbitrary real normalization constant A drops out), 
we obtain e, and therefore the critical temperature 

^ = 4;^12V4 ' 

as a function of the charge q. The solution may be written in terms of the hypcrgeometric functions in (|35p . but we 
will not present it here. At the minimum value = qc — found from the critical line ^ with = —2, we 
obtain Tc = 0. For q^ > q^. , Tc is an increasing function of q^. This behavior is shown in figure [T] 
For q^ = 0, we deduce from ([35]), 



e « 0.004 , (37) 

and the critical temperature is 

— wl.TxlO""^. (38) 

For the other boundary condition, A = 1, we have ^ ^ z + 0{z^) as z 0. The calculation may be done in a similar 
manner, giving the value at q^ = 0, 

e « .553 , (39) 

and the critical temperature 

— w .024 . (40) 
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FIG. 1. The critical temperature Tc as a function of found by a numerical solution (dashed) and an analytic one based on 
hypergeometric functions (solid) for A = 1 (left panel) and A = 2 (right panel). Tc —> as —j. 



The critical temperature as a function of is shown in comparison with numerical results Q in figure [T] It should 
be pointed out that our analytic results are valid in the small regime where e is small, since we are performing a 
small e expansion. As increases, so does e, and so corrections from higher perturbative orders become important. 
Also, comparing different scaling dimensions, the corrections for A = 1 (eq. p9p ) arc larger, because e is not as small 
as for (eq. p7p ). They become increasingly accurate as we approach the critical point = ^i- this limit, 

the numerical analysis becomes unstable. 



B. Variational method 



Next, we solve the wave equation and determine the critical temperature for arbitrary scaling dimension A using a 
variational method. Our results extend earlier work by HartnoU, et al. Q who used a variational method to show that 
there was a phase transition in the case = — 2, g-^ = 0. Their trial functions worked well for A = A_, however they 
were unable to find a trial function which would show instability in the case A = A_|_ where the critical temperature 
is lower. We find trial functions for all scaling dimensions A ~ A± and coupling constants that show that we have 
instability down to a critical line = ?c(^) that we compute analytically. 

Let 



(Oa> 
V2 ■ 



The equation of motion for the scalar field (|17l) in the background (|T5|) reads 



F" 



2(A-1) 



F' 



A 

72 



(A -3) 



.9 



F 



2 

q p — 



-F . 



For a given p, q^ is an eigenvalue determining Tc- It extremizes the expression 



9 = ^ 



V 



where 



A/" 



dz z 



2A-2 



\gF'^ + Az(A ■ 



A^(l 

4 ^ 



-z)F'^\ , V 



dz z 
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(41) 



(42) 



(43) 



(44) 



Of the (generally infinite number) of extrema (eigenvalues), we need to select the minimum. To this end, consider the 
one-parameter family of trial functions 



(45) 
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where Z- is the inner horizon (z^ > = 1 and g{z-) — 0). 

The critical point can be found by substituting the trial function (|45|) in (|44|) . and taking the extremal limit 
12. We obtain 

3-2A(3- A) 



JV = ^ ln(z_ - 1) + na{a) + Oi^J~T) 



V = -2o? ln(z_ - 1) + (io(a) + z^ - 1) . (46) 

The functions no (a) and do (a) can be found explicitly, but will not be needed for our purposes. 
In the extremal limit z_ — > 1, we obtain from eq. (|43p . 



, _ 3 - 2A(3 - A) , / 1 



and therefore g^, which is given by (j47p in the limit z„ — > 1, is in agreement with expectations from geometrical 
considerations at the horizon (eq. ([1])). 



The above conclusion is valid as long as a ^ 0. In general, setting a = 0, one obtains a value of the ratio (|43|) 
which is higher than the minimum, validating the above conclusion. However, for sufficiently small A, the value of 
the ratio (|43p at a = is lower than the value (l47l) . In this case, the minimum is attained at a = 0. We deduce the 
critical point (in the limit Z- — >■ 1), 

2 ^ no(0) ^ A(1~A)(2A-1) ^^o , . 

do(0) 4(2A + 1) 3i^F(2A-l,l;2A;^) ■ ^ ^ 

The two possible critical coupling constants are plotted in figure [5] as functions of the scaling dimension A. The 
two critical lines meet at Ag which solves 

3 - 2Ao(3 - Ao) _ Ao(l - Ao)(2Ao - 1) Q^o ^^^^ 



Numerically, 



4(2Ao + l) 3^f^(2Ao-l,l;2Ao;^) 



Ao « 0.64 . (50) 



For A > Ao, the critical point is given by ([T]), whereas for A < Ao, it is given by (|48p . We have instability in the 
shaded region of figure [2] (above the minimum of the two critical lines for each A). Notice that even in the range 
A < Aq, we have < 0, with = at the unitarity bound A = 0.5. Therefore Reissner-Nordstrom black holes 
are unstable against neutral hair down to the unitarity bound, contrary to what one would expect by geometrical 
considerations at the horizon. This result is confirmed by a numerical calculation of the critical temperature as a 
function of for various values of A (figure [S]) . 

The above results show conclusively that we have instability for all > q^. Indeed, for a given q^ > q^^ define the 
action 

S=N-q^V, (51) 

with M and T) given in (|44p . The action S is extremized {6S = 0) by the solutions of the wave equation (H^ . It 
becomes negative for trial functions corresponding to eigenvalues (coupling constants) < q^, since 

S^{f - q^)V < . (52) 

This shows that we have instability above the critical point {q^ > Qc)- -By the same argument, below the critical point 
{q^ < q^), we have S > 0, since if > q^ > q^ for all trial functions, therefore no instability against scalar hair. 

We obtain analytic expressions for the critical temperature as a function of q^ for given A by using These 
expressions are compared with numerical data in figure |4l We obtain excellent agreement (the curves are almost 
indistinguishable) both near and away from the critical point, as shown in figure |H The variational method is 
superior to the numerical solution as we approach the critical point (Tc — >■ 0), because the latter becomes unstable 
there. 



FIG. 2. The critical coupling constant vs. A (boundary of the shaded region). The curve close to the horizontal axis 
represents the result H48|) of our variational method, whereas the steeper curve is the result ([T]) one obtains from geometrical 
considerations. 



IV. BELOW THE CRITICAL TEMPERATURE 



Below the critical temperature ([22|) . ^' does not vanish. We shall use perturbation theory to obtain solutions to 
the equations of motion. We expect the pcrturbative expansion to converge rapidly because the critical temperature 
is very small near extremality. We therefore expand 

(Oa)*o + (Oa)^«'i + ... 
$ = $0 + (Ca)^*! + • • ■ 
X = Xo + (Oa>'xi + ... 

9^90 + {Oa?9i + ■■■ (53) 
The charge density (in units in which the radius of the horizon is = 1) is similarly expanded, 

P=-$'(0) = po + (Oa>V + --- (54) 

where po is the value of the charge density g in units of the radius of the horizon at the critical temperature. 
The temperature of the hairy black hole (|13p is expanded as 

^ = ^[1-(Oa)Vi + ...] , (55) 



around the critical temperature (|22p . where 



At zeroth order, we have 

$o = Po(l-^), Xo = 0, g, = l-(^l + ^jz' + ^z\ (57) 
and 5*0 obeys the wave equation in the background ([57)1 and has already been calculated, 

^0 = ^^^^(^) ' (58) 
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where F is normalized by F{Q) = 1. 
At first order, we obtain 



Xi{z) = / dz' 
Jo 



We deduce the correction to the electric potential 

^ - r dz'^[{z') , $;(z) = -Pi + po^iz) 

J z 

The parameter pi is an integration constant to be determined, and 



ip{2 



dz' 



l-z' ^ z'(\-z'f 



/\\2 



(*o(/)) 



We also obtain 



where ^(1) =0 and 
Therefore at the horizon. 



2z4 2z2go(z) 



+ ^(*oWf + f (Xi W - 2^(z)) 



(59) 



(60) 



(61) 



(62) 



(63) 



(64) 




-0.010 -0.008 -0.006 -0.004 -0.002 



FIG. 3. The critical temperature vs. the coupling constant for A — 0.634, 0.578, 0.525, 0.501 (curves become steeper as A 
decreases) . 
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FIG. 4. The critical temperature Tc found by a variational method (dashed) compared to numerical results (solid), for A = 1 
(left panel) and A = 2 (right panel). 



The remaining unknown parameter pi is found from the first-order correction to the wave equation, 



"50 


2 




.50 


z 





.9o ^^90 



where 



9'\ gig'o x'l 



_9o 9o 2 

After taking the inner product with ^f'o and using the zeroth-order wave equation, we arrive at the condition 

1 dz 



determining pi. This is linear in p\ and easily solved, 



Pi 



Jo 1 (2 c^j 50 — 



^zG , 2q^ pl(\-z)(go^-z\\-z)g) 



90 



*0 



1 



^^^(3-42) _ g'oz3{l~z) 
2 2go 



*o«'o + 



rn^z{l-z) _ q-^il-zyCIgo+plz^l-z}} 
2so 



*2 



(65) 



(66) 



(67) 



(68) 



From the second-order expression (j55p for the temperature, we deduce the leading behavior of the condensate below 
the critical temperature. 



By setting T = 0, we deduce the first-order approximation to the energy gap, 

(Oa)1/^ 1 



1/(2A) 



(69) 



(70) 



In figure [5] we compare our first-order analytic results for the gap at = to numerical results obtained by solving 
the equations of motion numerically. The agreement improves as one decreases q^, because the critical temperature 
Tc decreases. In figure [51 we plot the gap as a function of and show that it diverges as — )■ (i.e., Tc — > 0). 
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FIG. 5. The energy gap for — 0, and A = 1, 1.5, 1.7 (left to right). The thin hne is our first-order analytic result, the thick 
one is from numerics. 
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FIG. 6. Left panel: energy gap vs temperature for A = 1.5 — 0.38 (lower curve), q^ = (middle curve), and q^ = —.125 
(upper curve). Right panel: energy gap vs q^ at temperature T = .95Tc (asymptote at q^ = —0.375). 



V. CONDUCTIVITY 



Next, we calculate the conductivity of the system. To this end, we introduce perturbations for the electromagnetic 
gauge field while respecting translational invariance via 

A, ^ gt. ^ e-*"* . (71) 

The equations of motion are Q 

K9e~) +{— 2q^—e—\A, = -^'e-{z^g't^ + 2zgt.,) , 

2g,, + z{<i>'A,,+g',J=0 , (72) 
which yield the equation for the electromagnetic perturbation (in Schrodingcr form) 

- ^ + = ' (73) 

dzi 

where is a "tortoise" coordinate defined by 

"'.eX(^')/2 ^ (74) 



" dz' 



5(^0 



with = at z = 0, and the potential is given by 



V ^ g i 2q'''^e-^ + z^<^>'^ ] . (75) 



z 
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The wave equation ([73]) ought to be solved subject to ingoing boundary condition. ~ e'"^* , at the horizon. At the 
boundary, expanding the gauge field in a series, 



A, 



we read off the conductivity, 



4°' + zAi'^ + 



iujA 



(0) 



(76) 



(77) 



To visualize this, it is convenient to extend the domain to the entire real axis with F = for < 0. Consider a wave 

,iw2, Jqj. ^ Q_ p^^ _ gg^g partly reflected and partly transmitted. 



of "energy" w incident from the left, A^ 



so that for < 0, A^ 



whereas for > 0, we have a transmitted wave which is purely ingoing at 



(78) 



the horizon. The reflection coefficient is then related to the conductivity (|77p by jlc 

1 -7^ 



1 + 7^ 



To find IZ. write the wave equation in integral form as 



dz:G(z,,z:)y(z:)A,(z:) , 



where the Green function satisfies 



Explicitly, 



G(z„z:) = — 



The reflection coefficient is then 



and can be calculated perturbatively, 



7^ 



z,)e*"("--"*) + 0(z* - 2;)e-*"«-^*) 



dz*G'(0, z^)V{z^)A^{z^^ 



7^ 



dz,G(0,zO^(^*)e*'^^- + / dz,G(0,z,)V(z,) / dz',G{z,,z',)V{z',)e^'^''' + 



(79) 



(80) 



(81) 



(82) 



(83) 



The numerical results are shown in figure [T] The analytic results at second order are compared with the numerical 
solution in figure[51 As expected, we have agreement at high cj. As we go to higher perturbative orders, the agreement 
extends to a wider range of frequencies. 

We can look at the potential to determine the behavior of the conductivity near the critical line ([T]). We shall 
concentrate on the case = 0, in which the system simplifies considerably. Our discussion can be extended to other 
values of the coupling constant in a straightforward, albeit tedious, manner, but no new physical insight is gained. 

Regardless of the value of A, the potential reaches a universal shape and never becomes large enough for spikes 
to form at or near the BF bound, as seen previously in the probe limit [l^, [13, [13 ■ The potential is plotted in 
figure ini below the critical temperature. After scaling V by the charge density q and z — >■ Zy^, the curves at various 
temperatures T < Tc are almost indistinguishable. Detail near the horizon (right panel of figure [H]) shows the effect of 
lowering the temperature: the potential approaches zero faster at the horizon resulting in a longer tail. Consequently, 
the conductivity also approaches a universal shape (figure [7]) below the critical temperature. 



VI. CONCLUSION 



We discussed the properties of a hairy black hole near extremality corresponding (via a gauge/gravity duality) to 
a system of low ratio Tc/ fi where Tc is the critical temperature at which an instability arises and /i is the chemical 
potential. Using hypergeometric fucntions in the special case = —2, and trial functions in general, we calculated 
Tc analytically for small charge of the scalar hair continuing analytically to negative values of as needed, down 
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Im(o-) 



-0.2 



5 E, 




5 E, 



FIG. 7. Real and imaginary parts of the conductivity. 




5 E, 




FIG. 8. Comparison of analytic and numerical results on the real and imaginary parts of the conductivity (left and right panel, 
respectively). 
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FIG. 9. Left panel: the potential determining the conductivity for q = at various temperatures below Tc. Right panel shows 
detail near the horizon. 



to Tc = 0, which occurs at a critical value = q^. Below q"^, there is no instability and no scalar hair develops. 
The critical point q^ was expected to be given by ([1} as a function of the scaling dimension A. which is of geometric 
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origin, arising from the AdS2 geometry of the horizon at extremahty. We showed that this is only true for A > Aq, 
where Aq ~ 0.64, whereas below this value, the true critical line is lower. As we approach the unitarity bound, the 
critical coupling remains negative, with = at A = 0.5. Among other things, this implies that for q ^ 0, we have 
instability all the way down to the unitarity bound contrary to geometrical expectations. 

Below the critical temperature, we calculated the condensate and the conductivity. The energy gap in units of Tc 
was foimd to diverge as we approached the critical line (i.e., Tc — ^ 0), and the conductivity as a function of frequency 
(appropriately normalized) was found to have a universal shape below the critical temperature. In particular, there 
is no strong dependence on the charge, which is different from the probe limit where there is a frequency gap that 
has an explicit q dependence {ujg ^ (qOA)^^^)- 

Despite the mild dependence of physical quantities on the charge, it should be noted that the system of zero charge 
has an enhanced symmetry, since the phase of the scalar field decouples. It is therefore worth exploring its properties 
further, perhaps using a solution to the equations of motion with a non-trivial phase dependence, such as a vortex. 
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